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CHAPTER I 


INTRODUCTION 

The prediction of temperatures and temperature gradients within a 
spacecraft while it is operating in a space environment has become a 
significant problem m the design of such spacecrafts. Performance tests 
on full sized prototypes of the spacecraft m a simulation chamber were 
first used to solve this problem but the advent of extremely large space- 
crafts and space stations has made this technique become outmoded. As a 
result, thermal modeling has reached a level of increasing importance, 
since it can successfully be used to predict the thermal behavior of a 
spacecraft from model test data. 

A great abundance of analytical and experimental work has been done 
in the area of thermal modeling for steady state conditions imposed on 
conduction-radiation coupled systems m recent years. The transient be- 
havior of thermal models has been greatly neglected, however, and the 
need for more work m this area is readily apparent. 

The primary objective of this investigation was to develop a practi- 
cal method for thermal modeling of a conduction-radiation design under 
transient conditions. The design selected for this investigation consisted 
of solid rods. The models and prototypes were tested in a cold, high- 
vacuum environment. Simulated solar or particle radiation was not attempted. 
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CHAPTER II 


LITERATURE SURVEY 

In 1962 the results of thermal modeling of spacecraft and thexr com- 
ponent parts began to appear in the literature, Vickers (1) published a 
literature review of thermal modeling efforts covering the years 1962-1965, 
Vickers emphasized that while many text and papers on thermal similitude 
and modeling had been written, most of these mainly concerned modeling of 
convective problems and payed little or no attention to the modeling of 
conduction or radiative situations. This is indeed unfortunate since coupled 
conduction and radiation predominates m spacecraft, 

Vickers discussed two basic methods of thermal modeling The first 
method was "temperature preservation", where the model and prototype abso- 
lute temperatures are equal. In the "material preservation" method the 
model and prototype are constructed of the same material, Vickers concluded 
that the two methods were mutually exclusive, assuming perfect geometric 
similarity between the model and prototype. The article also discussed 
general problems of thermal modeling, reasons for improving the state-of- 
the-art, and advantages 'to he derived therefrom. 

As Vickers stated, the first paper m the area of spacecraft thermal 
modeling was by S, Katzoof (2). Katzoof derived five dimensionless groups 
from dimensional analysis, and the major portion of his paper was devoted 
to the practical problems m using the groups he presents. The groups 
Katzoof derived were concerned with radiation, internal heat generation, 
thermal conductivities of materials, heat capacities of materials, and 
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joint conductance. The five parameters were: 

oT^ , » cjL , C p ^ , and CL (1) 

4 <1 kT 4 6 k 

The last ratio m equation (1) applies to joint conductance. This group 
of ratios does not explicitly include radiation exchange between surfaces ; 
however, they form an independent group. If joint conductance is neglected, 
equation (1) becomes: 


CT 


c LT 

_P » > and S_L 

q q e k T q 


( 2 ) 


In 1962, A. J. Katz (3) discussed very briefly the criteria involved 
m thermal scale modeling, and some of the problems which might be involved. 
Katz, Clark, and Leband (4) were the first investigators to use dimensional 
analysis for establishing the thermal modeling criteria m space applications, 
and a very detailed paper was presented by this group. They derived a group 
of dimensionless ratios used to study the temperature of composite wall 
constructions exposed to an external radiation environment. Leband and 


Clark presented the ratios m equation (3) 

.2 - 


T c L 
dL p 


aT 


q ’ kT ’ k9 


(3) 


and Katz presented the ratios m equation (4) 
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3 

oT L 
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(4) 
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The two groups of ratios m equations (3) and (4) may be combined to form 
equation (5) : 


cjL 2 oT 3 L C p L 

k T ’ kT » k ’ k6 



(51 


The ratios in equation (5) that form an independent group are 
2 


3 

oT L 


c L 
P 


4 L 


, and 


q L 


k 5 k 0 * k T ’ “““ k T 

If joint conductance were neglected, equation (1) would become 


£l- 3k and 3Jt 

q * q 6 * kT » 4 


( 6 ) 


(7) 


Equations (6) and (7) are not identical groups of independent ratios, 
although they apply to the same physical problem. Either of the groups 
may be derived from the other. 

Wamwright, Kelly, and Keesee (5) discussed the development of dimen- 
sionless groups from both a dimensional analysis approach and from a dif- 
ferential equation approach. The similarity groups obtained were as 
follows : 


k 

o 
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c L 
o po 
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aeL T 
o 

k (AR) 
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( 8 ) 


Hrycak and Unger (6) developed the same thin shell thermal modeling criteria 

that are presented m equation (8) when the system is referred to m spherical 
coordinates. Such prototypes as Telstar were modeled by these investigations. 
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In dimensional analysis the pertinent variables are identified at 
the outset. Some insight into the problem and knowledge of at least one 
combination of parameters are considered essential to the success of this 
method. The entire group of dimensionless combinations of these para- 
meters formulates the modeling criteria. 

The similitude approach consists of writing the governing differential 
equation for the model and prototype thermal behavior. The model equation 
is identical to the prototype equation if the dimensional elements multi- 
plying corresponding terms are identical m value. The modeling criteria 
are established by equating the ratios of corresponding dimensional multi- 
plying elements m the model and the prototype differential equations. 

B. P, Jones (7) and (8) used the differential equations which des- 
cribed the thermal behavior of a spacecraft to derive the similarity para- 
meters of thermal modeling. Jones developed twenty-eight dimensionless 
ratios that had to be constant from prototype to model if complete simi- 
larity was to exist. All of the ratios were not independent, but it was 
possible to find seven independent groups of ratios. A comparison of his 
resulting similarity parameters with those of Vickers and other investi- 
gators was made by Jones and the results were found to be the same. 

Chao and Wedekind (9) assumed a relationship between thermal con- 
ductivity and temperature, ignoring only polarization and assuming in- 
coherency of the light, and m their approach state the necessary and 

sufficient conditions for modeling real surfaces under any form of solar 
input or simulation. They introduced the concept of using a pox^er lax* 
function of temperature for specific heat and thermal conductivity and 
also discussed the possibility of the model not being geometrically similar 
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to the prototype. 

Fowle , Gabron, and Vickers (10) made an elaborate attempt at scale 
modeling m order to predict the temperature of the Mariner IV spacecraft. 

The Mariner IV was exposed only to external radiation from the sun, and 
was at a fixed altitude with respect to the sun. This allowed steady 
state modeling to be applied. A prototype, which simulated the basic 
features of the Mariner spacecraft, was built and was then modeled at 
one-half and one-fifth scales using the the temperature-preservation 
technique and at one-half scale using the materials preservation techi- 
que. Reasonable success was obtained with all three models; that is, 
discrepancies between model and prototype were, m almost all cases, 
accounted for. The models were tested at steady state with cold-wall 
and high-vacuum conditions. 

In 1963 Vickers (11) again discussed the techniques which could be 
evolved from the basic laws of thermal scale modeling of spacecraft. The 
equations presented for the similarity parameters by Vickers m this article 
agreed with those developed by Fowle, Gabron, and Vickers, Jones, and earlier 
writers. The dimensionless groups were. 


nT 2 pc L 
qL P 

kT ’ k0 


m 3 T a, <f» i h a SL 

eaT L Cl kj ik s 

’ k " * k 5 kT ’ kT ’ 


Sk 

LkT 5 kT 


(9) 


If q, q, and q are considered as distinct physical quantities, equation 

(9) forms one group of independent ratios. The first ratio and last two 

ratios listed m equation (9) may be reduced to one by introducing the re- 
2 3 

lations q = q/L and q = q/L . The resulting group contains six independent 


ratios which includes one for joint interface conductance. In this case q is 
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a characteristic power; however, in certain modeling it may be desirable 
to consider q, q, and q as distinct quantities. 

In 1964 a paper by J. R. Watkins (12) revealed that he had 
developed forty-nine groups, each of which contained seven independent 
similarity ratios for the general case of thermal modeling. Watkins 
used the dimensional analysis approach to develop similarity ratios. 

By defining the energ y transfer to and from single, elemental, isothermal 
volumes of the prototype and model m a simulated space environment, 

Watkins was able to use the computer to derive his similarity ratios. 

A numerical approach to dimensional analysis was also applied to the 
physical quantities that describe the thermal behavior of two radiating 
disks connected by a conducting rod* The numerical solution yielded fifty- 
seven groups each of which contained five independent ratios. Any one 
of these groups may be used for model design, the selection depending 
upon what is to be the purpose of the model experiments. 

Matheny (13) was one of the earlier investigators to conduct experi- 
ments that involved transient thermal modeling. For his system Matheny 
chose two disks connected by a conducting member. The conducting member 
was made large enough to enter into radiant interchange. The models were 
designed according to the criteria 

•A c L 2 . . t 2 

al L p qL and qL 

k ’ k0 ’ kT ’ kT ^ } 

The results of the transient tests which were conducted on the prototype 
and one-half scale model were excellent. 
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The modeling criteria given by equation (8) were applied to a con- 
ceptual space station by Folkman, Baldwin, and Wamwnght (14) . Geometric 
scaling was used and external radiation sources were preserved. The ex- 
perimental results obtained from tests performed on the model were com- 
pared with three-dimensional transient analysis. 

R. E. Rolling (15) and (16) developed modeling criteria for space 
vehicles by using the similitude approach. Rolling's procedure was identical 
to that of Chao and Wedekind, The criteria m equation (11) were developed 

for prediction of prototype behavior from model behavior for a radiation- 
conduction coupled problem, where the "starred 11 terms represent the ratio 
of model to prototype properties. 


p*V*c T* 
P 

6 * 


k k A * k * 

A S* F - A R* F = A E* F 
s s r r e e 

* * * * 4 

q- = k A (T-) = A F T* 
n n X* J J 


( 11 ) 


The criteria in equation (11) were used to model two opposed disk with 
four connecting tubular members and truncated cones. 

Rolling used one-half and one-quarter scale models m his tests, and 
fixed material thermophysical properties while temperature and time were 
scaled. Geometric distortion in the minor dimensions was permitted. 
Arrays of tungsten filament lamps were used for external radiation 
sources. The results for the opposed connected disks and for the cones 

proved that the criteria m equation (11) were reasonably accurate for 
modeling purposes. The temperatures of the model for the disks were 
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within 9°K of the prototype, and for the cones all model temperatures were 
within 8°K of the prototype. 

Young and Shanklin (17) modeled a system composed of three sections* 
a heater plate containing an electrical resistance heating element which 
represented an internal energy source of variable strength, a cold-plate 
having no internal energy source, and a cylindrical rod connecting the 
two plates. The ratios m equation (12) were applied to this system. 




» R R R 
„2 _r_ r 

q ’ R » L L 


„2 » R 
R r 


and 


l 2 f 




R 


( 12 ) 


The prototype and the one-half scale model were fabricated of the same 
material with identical values of e and a. The prototype and the model 
were exposed to identical simulated space conditions and were started at 
the same initial temperature, T . Furthermore, the prototype and model 
surfaces were blackened such that e and a approached unity. The ratios 
in equation (12) were used m designing the model. The system used and 
a large portion of the results of the investigation are shown in Figure 
1 and Figure 2, respectively. Model temperatures deviated from those of 
the prototype by an average of approximately 1.5 percent. 

Gabron and Johnson (18) submitted a report to the Jet Propulsion 
Laboratory describing their work on thermal modeling of the Mariner Mars 
64 spacecraft. A forty- three hundredth scale model of the spacecraft was 
fabricated and tested under steady-state thermal conditions. The "temp- 
erature preservation" technique was employed, with no solar radiation. 
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MODEL GEOMETRY PROTOTYPE GEOMETRY 


FIGURE I SPOOL CONFIGURATION 


{SOURCE’ REF 17) 




FIGURE 2 TEMPERATURE- TIME CHARACTERISTICS OF THE SYSTEM 

SHOWN IN FIGURE I (SOURCE: REF. IT) 
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It was concluded that "temperatures within a complex spacecraft structure, 
typical of the Mariner vehicle, can be predicted by use of 'temperature 
preservation' thermal-scale modeling techniques to an accuracy useful for 
the verification of thermal design concepts." 

Jones and Harrison (19) used the group of ratios given m equations 
(13) and (14) to model a system composed of a plate, cylinder, and sphere 
which were exchanging thermal energy by radiation only. 



£ k,i F ki A ^c 
E A. 


a S A S6 
- 2 — 3 

c T 
3 3 


3 

E oA T 9 
3 3 3 

c 

3 


c, e 
k 3 


q_ 6 

and - — J 


c T. 
3 3 


(13) 


(14) 


Electrical resistance heaters were used to obtain the simulated heating 
effects of the sun, and the space chamber was regarded as incorporated 
m the ratio involving radiative interchange between surfaces. Equations 
(13) and (14) then reduce to 


T 

JL 

T,_ 


K T 3 

KJ 3 
c 

J 



% 9 

— 1 (j, k = 1. ..., n k t - j) 

c T 3 

3 3 


(15) 


The materials were not changed from prototype to model, the radiation 
geometry factors were assumed to remain unchanged, the temperatures of the 
model at a particular time were assumed to be equal to the corresponding 
prototype temperatures at the same time, and the thickness of the cylinder 
end caps was not changed from prototype to model* The major external di- 
mensions were scaled by one-half. The experimental results generally 
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confirmed the modeling rules. 

Adkins (20) introduced a method of geometric scaling which made 
thermal modeling possible while preserving both material and temperature. 
The Fourier conduction equation with radiation boundary conditions was 
used to establish the similarity criteria for two different configurations, 
a nonconduction array and a conduction array. The similarity ratios were 


T 2 

L pc 
P 

ke 


L 

d 5 


(crj d)/k 


(16) 


In 1966, Miller (21) investigated the applicability of thermal 
modeling to steady-state and transient conduction m cylindrical solid 
rods, for both single and multiple material systems. By performing a 
similitude analysis, Miller obtained the following relations: 

q* = L* -2 , q* = L* _1 , q* = L* 3 , T* = 1, 0* = L* 2 , R* = L* 2 

where the starred quantities again represent the model to prototype 
ratio of the parameter. 

Shih (22) presented a discussion of thermal modeling which included 

the presentation of some of the same similianty parameters as Miller. He 

took into account changes of thermal conductivity and specific heat with 

2 

temperature and included the geometric relationship R* = L~ for preserving 
both material and temperature m model and prototype. No experimental work 
was presented or planned, nor was there any reference to joined materials 
with different thermal properties. 

Price (23) investigated a model that maintained both material and 
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temperature preservation at homologous locations. The prototype and 
model, which were built and tested, were cylindrical fins. Also, math- 
ematical models of the prototype and model were solved using the computer. 
The similarity criteria used were those developed by Miller. Experimental 
results showed that the scaling laws were valid. Also the theoretical 
data agreed with enough accuracy to validate using the mathematical 
model, when the model is simple. 
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Summary 


Analytical and experimental research m thermal similitude m 
space-related thermophysics problems was very active from 1962 to 1966, 
hut has become somewhat inactive since 1966. However * a rejuvenation 
of this research occurred m 1968, and more emphasis is being placed on 
thermal modeling endeavors. Much progress has been made m confirming 
the analyses by experiments m both steady-state and transient cases. 
However, more work is needed m the conduction-radiation area with em- 
phasis on transient problems. There is also a great need for verifica- 
tion of a group of similarity parameters for a system containing a fluid. 
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CHAPTER III 


MODELING CRITERIA 

The obvious objective of thermal scale modeling is to maintain ther- 
mal similitude between a prototype and model. This is achieved by pre- 
serving the numerical equality of certain dimensionless groups of pro- 
perties of the prototype and model. These thermal similarity parameters 
may be identified from either dimensional analysis or from dependence 
upon the differential equations describing the thermal behavior of the 
model and/or prototype. Eor the inexperienced, the latter method of 
using the differential equations and boundary conditions is much more 
desirable. However, either method will result m the same set of similar- 
ity parameters. 

The modeling criteria for the models and prototypes considered in 
this investigation were developed using the differential method. This 
technique was found m the reference by Murnll, Pike, and Smith (25) 
and was employed m getting the similarity parameters. Before proceeding 
to the derivations of the parameters, some discussion is m order con- 
cerning the constraints placed on the problem to facilitate the derivation 
process. 

The first restriction was that the materials used be homogeneous 
and isotropic. The second basic assumption was that there be perfect 
geometric similarity between prototype and model. Third, to eliminate 
the influence of spectral and angular distribution on the emitted radia- 
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tion, it was assumed that the prototype and model had the same uniform 
and constant surface characteristics. To assure control of the third 
constraint, the surface of the prototypes and models was coated by a 
flat black velvet paint, which in turn assured the emissivity and 
absorptivity of the surfaces be constant and independent of temperature. 

It was also assumed that the radiant heat flux from a simulated space 
environment is approximately zero, which is a safe assumption if solar 
energy is not simulated and the temperature approaches absolute zero. 

As a consequence of the fourth assumption, it was also assumed that all 
energy radiated from the models and prototypes was absorbed by the walls 
of the cryogenic liner inside the simulation chamber. A fifth assumption 
was that thermal properties of the model and prototype were constant and 
did not vary during a test run. To assure this, a temperature range of 
approximately thirty to ninety degrees Fahrenheit was used during all of 
the test runs. Finally, heat transfer due to convection was assumed to 
be negligible, which is an excellent assumption m a vacuum environment. 

With these constraints in mind, the decision was made to test titan- 
ium and stainless steel rods to verify the modeling criteria for a con- 
duction-radiation problem. These rods were very long relative to diameter 
for the purpose of minimizing end effects. The technique found in re- 
ference (25) was used to develop the similarity parameters. 

In order to derive the similarity parameters needed for thermal scale 
modeling, it is'" necessary to make appropriate changes of variables m the 
describing differential equation. A step-by-step technique was presented 
by Murrill, Pike and Smith for making the correct changes of variables 
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for two basic reasons * First, those inexperienced m the field could 
use this technique to derive the proper changes of variables* Second, 
this systematic technique could be used for problems m which the proper 
changes of variables are not .obvious* Normally, the changes of variables 
are made for the dependent and independent variables m order to reduce 
the number of groups of parameters required to adequately describe the 
problem. However, m a few cases the changes of variables are even more 
helpful. Sometimes the number of independent variables required m the 
equations can be reduced, such changes m the independent variables being 
known as similarity transforms* The derivation of similarity transforms 
is a degree more difficult than the derivation of changes of variables, 
but the technique presented will obtain such similarity transforms, if 
they exist, m addition to changes of variables* 

The techniques for reducing the number of parameters or dimensionless 
groups required to mathematically describe a physical problem have been 
investigated for many years* The description of problems using the mini- 
mum number of parameters is the subject of such apparently diversified 
fields as dimensional analysis, inspectional analysis, and scale-up theory. 
However, these techniques typically consider only changes of variable, 
and not similarity transforms* Since the following technique applies to 
similarity transforms as well, it would seem to be preferable* 

The basic steps m this method of analysis are as follows: 

1* Introduce arbitrary reference quantities m order to place 
the dependent and independent variables m dimensionless form. 

This typically takes the following form: 
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variable-arbitrary reference value 
arbitrary reference difference 

For example , the temperature T may be non- dimens lonaly zed as 
(T-T*)/AT*, where T* is some arbitrary reference value and AT* 
is some arbitrary characteristic temperature difference m the 
problem. 

2. Introduce these changes of variable into the differential 
equations and their boundary conditions, rearranging as necess- 
ary so that the parameters appear in dimensionless groups. 

3. Equate each dimensionless parameter containing an arbitrary 
reference quantity to either unity or zero until all arbitrary 
constants have been defined. 

In the third step a system of algebraic equations is obtained that 
relate the arbitrary constants to the original parameters of the problem. 
There are three possibilities for this set of algebraic equations. 

1. More equations than unknowns, i.e., an over determined set of 
equations. In this case there will be groups of parameters or 
dimensionless groups remaining m the transformed equations. 

2. Number of equations equals number of unknowns. For such situa- 
tions no dimensionless groups of parameters usually remain in 
the transformed equations. 

3. Fewer equations than unknowns, i.e*, an undetermined set. In 
this case, a similarity transform is possible. 

Derivation of Similarity Parameters 

The general heat-conduction equation governing the temperature 
distribution and the conduction heat flow in a solid cylindrical system 
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having uniform physical properties is given by Kreith (26) as, 

3 2 T 1 3T q gen _ pc 3T (18) 

. 2 r 3r k k 3t 

dr 

where q represents internal generation due to resistance heating, 
gen 

2 

and is equal to I (Res) * The initial and boundary conditions for the 

V 

rod m a low temperature vacuum environment are as follows: 


T(r,o) = T 

0 

(19) 

3T(o,t) n 

3r 

(20) 

-k = oe [t (R, t)] 4 

(21) 

The following dimensionless reference quantities will 

be introduced 

into equation (18) and the initial and boundary condition 

equations . 

T " T rf 

AT _ 
rf 

(22) 

r-r * 
, rf 

r* ss 

Ar rf 

(23) 

t * = t_trf 
it , 
rf 

(24) 

Equations (22), (23), and (24) may be rearranged to give. 


T = T*AT . + T _ 
rf rf 

(25) 

r = r*Ar - + r , 
rf rf 

(26) 
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(27) 


t = t*At ,, 
rf 


+ t 


rf 


The partial derivative of each of the above equations may next be found. 


AT 

3T = rf 3T* 

3r Ar _ 3r* 

rf 

= AT rf 3 2 T* 
3 r 2 ( Ar rf > 2 9 ~* 2 


(28) 


(29) 


AT 

3T = rf 3T* (30) 

3t At . 3t* 
rf 


Equations (28) , (29) , and (30) may be substituted into equation (18) 
and into the initial and boundary condition equations. Substitution into 
equation (18) gives: 


_1 \ 3T* q gen (4r rP _ pc 8T* (31) 

r r£ 3r* k AI rf -k 4t rf »t* 

Substitution into the initial and boundary condition equations (19), (20), 
and (21) gives the following results* 



T* (r*. 


At 


, T - T , 
rf ^ _ o rf 


rf 


AT 


rf 


(31a) 


3T* 


( ~ r rf 

Ar rf’ 



3r* 


= 0 


(32) 
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3T* ( R-r rf ) 

irr » * 

rf 

3r> 


o E (Ar rf ) (AT rf )' 


/R-r _ \ . T _ 

■ , t*' + rf 


1 4 


Ar 


rf 


AT 


rf 


( 33 ) 


The dimensionless groups appearing m equations (31) through (33) are 
as follows: 


'rf 


(34) 


Ar 


rf 


q (Ar _)' 
^gen rf 


(35) 


k AT 


rf 


(Ar _)' 
pc rf 

k ( At rf > 


(36) 


rf 
At , 


(37) 


T - T - 
o rf 

AT - 
rf 

R-r- 

rf 


Ar 


rf 


(38) 


(39) 


rf 


AT 


rf 


(40) 


ae ^ Ar r f ^ ^ AT rf ^ 


( 41 ) 
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Since there are eight dimensionless groups m equations (34) 
through _ (41) , eight algebraic equations could be obtained of all these 
groups were set equal to zero or unity* However, this would be an over- 
determined set* There are six arbitrary constants; so, only six of the 
groups can be selected. 

The process of selecting the six groups is not as arbitrary as it 
might seem. The following reasoning was used to obtain the changes 
of variable 

1. The term t^ appears only m equation (37) and was set equal to 
zero to get the simplest form of equation (24) . 

2* The term T ^ appears m equations (38) and (40), and setting 

either of these equations equal to unity would complicate equation 

(22). Equation (40) was set equal to zero so that = 0* 

Setting equation (38) equal to zero would have given - T q , 

and this too is suitable since T was the same for model and 

o 

prototype. 

3. The r ^ term appears m equations (34) and (39). - Equation (34) 
was set equal to zero, which made = 0. This gave the simplest 
form of equation^ (23). 

4. Equation (39) was set equal to zero, which gave Ar^ f - R. This 

is the simplest form that can be derived for Ar^, and quarantees 
the presence of the rod radius m terms of R m the similarity 
parameters. 

5. Equation (36) is the only group in which At^ appears; so, this 
group was set equal to unity to define At 
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6. Equation (35) was set equal to unity to define AT^ and to make 
sure that the generation term appeared m the parameters. Power 
predictions for a prototype based on experimental power measure- 
ments from a model were used to achieve homologous temperatures 
between models and prototypes; so, it was indeed necessary to 
have a term for power (and m turn for generation) appear m the 
similarity parameters. 

This completed the selection of the six groups, and also the choice of 
whether to set each equal to zero or unity. 

At this point, then, the following expressions must hold: 
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At 


= 0 


rf 


(42) 
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AT 


- 0 
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(45) 


(46) 
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AT 


rf 
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Solving equations (42) through (47) for the six unknown reference 
quantities gives. 


t = 0 

rf 

(48) 

T _ = 0 
rf 

(49) 

r = 0 
rf 

(50) 

Ar = R 

rf 

(51) 

.,2 

At = 
rf k 

(52) 

q R 2 

AT = 

rf k 

(53) 


Substitution of equations (49) and (53) into equation (22) will satisfy 
the requirement of at homologous times. 




T k 
q 

gen 


R 2 


(54) 


For point to point similarity to exist between model and prototype, equation 
(55) must hold. 
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( 56 ) 



p m 

For equal thermal conduct xvx ties (model and prototype made of some materxal) 
2 2 

and q = (I Rcs/ttR L) - 
gen 9 


T 

m 

T 

P 



,I 2 Res 


P 


(57) 


Since T /T = 1 is 
m p 


required, equation (58) must hold: 


(tH 

\ 'm 


,I 2 Res/ 


P 
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(58) 


I = I 
p m 



Res /L 
m m 

Res /L 
P P 


(59) 


Equation (59) was convenient to use during the experiment to predict the 
current settings needed for prototypes. Equal starting temperatures and 
ending temperatures were achieved for models and prototypes using equation 
(59) . The approach, then was to check the scaling technique by using 
power predictions to achieve temperature equality. Of course, a time 
scaling factor was required m plotting temperature versus time curves 
for models and prototypes to prove that temperature equality had been 
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achieved at any instant m time. 


Substitution of equations (48) and (52) into equation (24) gave the 
time scaling equation between models and prototypes. 

t* - ~~2 (60) 

pcR 



Substitution of equations (50) and (51) into equation (23) gives: 



(63) 


or 


r 

in 

r 

P 


R 

P 


(64) 


Equation (64) has significance m that some value may be chosen for the 
radius ratio to be substituted into other modeling equations, i.e. 

R /R = -q-jr > etc. 

m p 10 

The following forms for the similarity parameters have been obtained: 
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(65) 
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( 66 ) 
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(67) 


2 

The (I Res) term is somewhat ambiguous for modeling purposes, but appears 
only because internal generation is a result of electrical resistance 
heating. Equation (65) may be put m a simpler form by substituting 
Res = p L/A, where P £ is the resistivity and depends on the material and 
temperature. There are numerous tables available that give the resistivity 
of common metals and alloys. Making the substitution for Res m equation 
(65) gives. 



( 68 ) 


Equation (68) assumes that the resistivity for model and prototype are 
equal and that = Res A/L accurately defines the resistivity. Tne ex- 
perimental program showed the actual resistivity to be approximately ten 
per cent less than the theoretical resistivity defined by the p^ equation. 

_This fact must be considered when equation (68) is used, plus the fact 
that if model and prototype are not fabricated from the same piece of 
stock the resistivities will differ slightly even if the model and pro- 
totype are made of the same material. 
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The final forms for the similarity parameters are as follows: 



A group of parameters were derived using equation (35), (36), and 
(41), since the radiation group was not included m the derivation above* 
The similarity parameters derived m this case gave reasonable results, 
but the results were not as accurate as those obtained by using equations 
(69), (70), and (71)* 

In most thermal modeling endeavors, the model to prototype ratio 
of the various parameters has a known constant value* This is parti- 
cularly true for the dimensions of radius, length 3 etc. Experimental 
knowledge of a model and intuition assist the investigator m choosing 
the value of this model to prototype ratio. 

For the rods tested m this investigation, the ratio (r^/r ) was 
known without making use of the scaling laws* Since radial heat conduc- 
tion and surface radiation were the only modes of heat transfer involved 
m testing the rods, the lengths of the rods were arbitrarily selected. 

As stated previously, the procedure was to test a model, and predict the 
power input for the prototype that would give (T^/T ) = 1 at homologous 
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points along the surface of the rods. If this predicted power input 
produced the desired temperatures, then the similarity parameters were 
known to be valid. The time scaling term m aquation. (70) could then be 
used to plot curves of temperature vs. time for the prototype and model 
to further validate the similarity parameters, as shown m Chapter VI. 
The curves should fall upon each other if the parameters are valid. 
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CHAPTER IV 


TEST MODELS 

The verification of the modeling criteria for a conduction-radiation 
model was accomplished by using metal rods, and this design was chosen 
for several reasons. The selection process was somewhat limited due to 
problems m having the test sections fabricated within a given time peri- 
od, The time factor also made it impossible to choose test sections that 
would lead to extremely lengthy theoretical work m developing the modeling 
criteria. Consequently, a solid rod was chosen as the most practical test 
design to check the validity of thermal modeling criteria. 

The solid rod met the requirements of giving one-dimensional (radial) 
heat conduction under transient conditions. The rod surfaces were finished 
at the same lathe speed to insure uniform and homologous surfaces on all 
test sections, which m turn insured uniform radiation from the rods. With 
the proper surface finishes it was then possible to achieve temperature con- 
sistency at each point where the thermocouples were attached to the rods. 

Two prototypes and two scaled models were constructed m order to 
check the validity of the modeling criteria. The models and prototypes 
were constructed from titanium alloy Ti-6AL~4V and type 304 stainless steel. 
Tables 1 and 2 give the actual dimensions of the prototypes and models used 
in the experimental phase of the program. 

Each of the rods was threaded approximately 0.50 inches on the ends 
so that the electrode terminal lugs could be firmly connected to the rods. 
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TABLE 1 


MODEL DIMENSIONS 


Model 

Length 
L (in) 

Diameter 
D (in) 

Titanium 

18.0 

0.255 

St. 304 

24.0 

0.2496 


TABLE 2 

PROTOTYPE DIMENSIONS 


Model 

Length 
L (in) 

Diameter 
D (m) 

Titanium 

24.0 

0.300 

St. 304 

24.0 

0.374 





















A "snug" connection was necessary to assure proper conduction of current 
through the rods. Other techniques of attaching the electrodes were tried, 
but these were found to be far too awkward and inefficient m that ex- 
treme heating often occurred at one end of a rod. 

As stated above, two electrical leads or electrodes were used to 
supply electrical power to the rods. The technique employed m sizing 
these leads is discussed in detail m Appendix A. One electrode was 
connected to each end of the rods. Two 20AWG copper leads (enamel coated) 
were also connected to each end of the rods to measure the electrical 
potential. The smaller leads were tightened around the threads by using 
pliers, and brass nuts were then screwed against the lead wires, which m 
turn pushed the wire firmly against the threads at the end of the rods. 
Following this the larger electrodes were attached to the rods. The ends 
of the electrodes were polished and were looped to fit the rod diameters. 
The electrode loop was slipped on the rod, and a nut screwed behind the 
loop forced it against the lead wire nut making an excellent electrical 
connection. 

The thermocouples used were 30AWG copper-const ant an, and the thermo- 
couple junction was made by soldering the wires together. This junction 
was then shaped to fit the curvature of the test section. The thermo- 
couple junction and lead wires were kept m place by using Devcon cement 
m the early stages of the experiment. However, this proved to be a poor 
technique and spot welding was used. Spot welding eliminated the possi- 
bility of the beads popping free of the rod during the test run, and 
minimized the response of the bead to voltage potential instead of 
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temperature at the various power settings. 

All thermocouple wires and the lead wires used to measure electrical 
potential were polished twelve inches from the point where they left the 
system. These wires were then extended and soldered to a junction plate 
kept inside the vacuum chamber. This junction plate was used to facilitate 
the changing of test models. The wires were extended from the junction 
plate to feed-throughs located m the chamber wall. All thermocouples 
were checked on a Leeds and Northrup potentiometer before a model or pro- 
totype was put into the 'chamber to prove them functional. 

Five thermocouples were attached to each of the models and prototypes. 
These were located at 16.66 per cent, 33.33 per cent, 50.0 per cent, 66.66 
per cent, and 83.33 per cent of the rod length, measured from one end of 
the rod. 

After the electrodes, electric potential leads, and thermocouples 
were attached to a rod, it was spray painted with a flat black paint, 
(Velvet coating 101-C10 by 3M) to a mean thickness of approximately 0.003 
inches. Nylon cord was then tied to the ends of the rods, and they were 
hung m the center of the vacuum chamber. 
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CHAPTER V 


EXPERIMENTAL PROGRAM 

In order to verify the modeling criteria developed m Chapter III, 
the models and prototypes described m Chapter IV were tested inside a 
space simulation chamber. The modeling criteria were used to predict 
the thermal behavior performed on the geometrically similar designs. The 
chamber provided the necessary low temperature-low pressure environment 
so that an accurate prediction of the energy exchange between the rods 
and the surroundings could be made. 

The space simulation chamber was a standard piece of test equipment 
located at NASA’s Mississippi Test Facility. The chamber was a Murphy- 
Miller high vacuum altitude test chamber made of stainless steel. It 
was forty-seven inches m diameter by sixty inches long less shelf height 
of four inches, and had one removable stainless steel flange. The test 
chamber was connected to a diffusion pump and roughing pump for producing 
high vacuums to one-half micron absolute pressure. The chamber was also 
equipped with feed-throughs for copper-constantan thermocouples, potential 
leads, liquid nitrogen, r and electric power. 

Since it was necessary to obtain a low temperature environment, a 
liner was placed m the mam chamber. The liner shell was constructed of 
stainless steel and was fifty-eight inches long and forty- three inches m 
diameter. The shell was spiral wrapped with five-eighth inch diameter 
copper tubing which carried the liquid nitrogen. At one end of the liner 
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was an optically tight baffle and at the other end was a removable flange. 
A liquid level controller was used m conjunction with a solenoidally 
operated valve to control the flow of nitrogen and to insure that liquid 
was m the tubing at all times during a test. The inner wall of the liner 
was coated with a flat black paint (velvet coating 101-C10) in order to 
insure surfaces with high and uniform values of emittance and absorptance. 
The outer wall of the liner and inner wall of the chamber were covered 
with foil to reduce the heat loss from the chamber. In order to minimize 
conduction from the outer chamber wall to the liner, it was necessary to 
support the liner on four adjustable legs. A schematic of the test faci- 
lity is shown m Figure 3, and a photograph of the overall vacuum system 
is shown in. Figure 4 (the fluid injection system shoxm m Figure 4 was 
used for experiments other than thermal modeling at MTF) . 

The inner wall of the liner reached -300°F, and the outer wall of the 
liner reached -275°F using liquid nitrogen. Figures 5 and 6 show the LN_ 
lines as well as hot water lines used to return the chamber temperature 
to room temperature as quickly as possible after a test run. 

To perform the transient tests on the rods, each rod was located cen- 
trally within the liner. The horizontal position of the rods was attained 
by extending nylon cords from hooks at the top of the liner to the ends 
of the rods. The nylon cord was very long m comparison to its diameter 
to minimize conduction losses. A photograph of the rod position is shown 
m Figure 7. 

The thermocouple output was obtained by using eight strip-chart re- 
corders. Five recorders provided the test rod temperatures, a recorder 
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gave the power lead temperature, another recorder provided the inside 
liner temperature, and the remaining recorder provided a measure of the 
current input to the rods. The recorders were calibrated with a Leeds 
and Northrup portable potentiometer with a reference temperature of 32°F. 
During the tests, the thermocouple reference function was maintained at 
32° F with an ice bath* Figure 8 shows the path followed by the thermo- 
couples m reference to the strip chart reocrders, i.e., the thermocouple 
circuit. Appendix B explains the procedure used to convert millivolt 
readings from the recorders to temperature readings. 

Before each test run the test rod was exposed to a vacuum environment 
at ambient temperature for a period of time. At this condition the out- 
put of all thermocouples was within + 0.20 millivolt of each other. It 
was therefore reasonable to assume that temperatues were measured to an 
accuracy of approximately + 1.00 during the test. 

For transient testing, a D-C power supply could be hand controlled 
to within + 1 millivolt during test runs. 

Current input to the rods was determined by measuring the millivolt 
drop across a calibrated shunt m series with the rods, using a vacuum- 
tube voltmeter. The VTVM was also used to measure the potential drop 
across the rods by making use of the leads attached to the rods for this 
purpose. 

The voltage drop and current input determined from the VTVM were used 
to check the resistance change of the rods as the power input was increased. 
The resistance was found to be relatively constant m the zero to one- 
hundred degrees Fahrenheit range used during the test. 
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Time measurement was accomplished by the time marker on the strip 

chart recorders. The error involved was negligible. 

The pressure inside the vacuum chamber was measured with a logatorr 

—6 

ion gage. The pressure was approximately 1 x 10 torr. 

After a test rod was installed within the liner, the thermocouples 
were connected to a junction plate, which m turn was connected to the 
power feed-throughs, the liner flange was installed, and the system was 
then ready for a test run. Pump -down was accomplished and the thermo- 
couples were checked for equality of output. When the chamber pressure 

^ 

reached approximately 1 x 10 torr, liquid nitrogen flow was started. 
After approximately three hours, complete cool-down was achieved. 

The copper tubing around the liner was maintained full of liquid 
nitrogen by a liquid level control system. The models and prototypes 
were protected from temperature extremes by adjusting the power input 
to the rods during the cool-down period. 

The input to the rods was set at a desired value and temperature 
recorded until steady conditions were established. Then the potential 
drop, current input, and temperature at the five thermocouple locations 
on the rod were recorded. A difference of approximately one-half degree 
per hour implied steady conditions. 

At this point, the power input was changed to a higher value and the 
temperatures and times were recorded until new steady conditions were 
reached. 
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CHAPTER VI 


EXPERIMENTAL RESULTS AND CONCLUSIONS 
Transient Results 

The experimental phase of the thermal scale modeling investigation 
was conducted using the facility described m Chapter V. The transient 
testing included testing of two prototypes and two models. The first 
test run was made using a titanium model and prototype, with R^/Rp = 0.825. 
Secondly, the stainless steel model and prototype were tested, with 
R^/Rp = 0.66. Each test was conducted m accordance with the procedure 
described m Chapter V. 

As stated m Chapter IV, the early part of the experiment was plagued 
by heat loss at the ends of the rods where the power leads were connected. 
Several test runs were made m an effort to eliminate this problem, and 
as stated m Chapter IV, a reliable technique was found. Several other 
miscellaneous problems were encountered in the earlier runs, but these 
were eventually eliminated. 

The energy input to the models was arbitrarily selected and a corres- 
pondingly scaled input was used for the prototypes according to the value 
predicted from equation (59). Temperatures were measured following the 
test procedures described m Chapter V. Tables III and IV show the scaled 
input values for the prototypes and the resulting temperatures obtained. 

Figures 9 and 10 present the measured temperature-time data for the 
two models and the corresponding time-scaled temperature-time data for 
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TABLE III 

TITANIUM RESULTS 
EXPERIMENTAL DATA 


FINAL TEMPERATURE 











TABLE IV 

STAINLESS STEEL RESULTS 
EXPERIMENTAL DATA 



INITIAL TEMPERATURE 

FINAL TEMPERATURE 


Model 

Prototype 

Model 

Prototype 

I (AMPS) 

15.30 

29.60 

18.65 

36.30 

E (VOLTS) 

, 0.236 

0.176 

0.280 


P (WATTS) 

3.611 

5.224 

5.222 

7.859 

TCI (°F) 

26.10 

27.50 

76.70 

78.00 

TC2 (°F) 

25.50 

27.00 

74.00 

77.50 

TC3 (°F) 

26.00 

26.50 

75.00 

76.00 

TC4 (°F) 

24.80 

26.00 

72.80 

75.40 

TC5 (°F) 

25.00 

25.00 

74.00 

76.40 
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the prototypes. Generally, there was excellent prediction of prototype 
temperatures from model data. The titanium results agreed perfectly, and 
the largest average absolute error for stainless steel model and proto- 
type was three and one-half degrees. A discussion of the errors involved 
is included m the following section of this chapter. 

Miller (21) conducted research with a solid cylindrical design made 
of two dissimilar materials. An interesting comparison exist between the 
similarity parameters Miller obtained and those obtained m this investi- 
gation. Some of the parameters he derived (neglecting property variations) 
were as follows: 




(72) 


T 

m 

T 

P 


1 


(73) 



(per) 
m 

<0 cr > m 


(74) 


These parameters were used for steady and transient predictions Since 

(pc) = (pc) equation (74) reduces to. 
m p 



The discovery was made that Miller's equation (72) predicted the same 
power settings as did equation (59), using the conversion 0.293 watts/btu/hr . 


50 



Consequently, temperature-time curves were drawn using Miller’s time- 


scaling equation (75). The results are shown m Fxugres 11 and 12, 

and indicate that use of equation (75) yields reasonably accurate results. 

Comparison of Figures (8), (9), (10), and (11) show that prototype temp- 

R 2 

ratures are slightly higher when time-scaled using t = t^ ^ as com~ 

R 

pared to Miller’s expression t - ^ . 1 

^ R 

m 


Error Analysis 

Some consideration was given to instrumentation errors m Chapter V. 
This section will discuss possible effects caused by such factors as in- 
correct power input to the model, energy losses through thermocouple leads 
and power leads, and variation of thermal properties. 

The power input to the models and prototypes was generally different 
than that required by the similarity parameters. Several test runs were 
made with each model and prototype, and the resulting quantity of data 

i 

provided a means of determining what effect slight power differences 
had on the rod temperature. The largest difference between actual power 
input and the input required by the parameters was +0.05 watts for the 
stainless steel prototype. This led to a +1 degree error, approximately. 
The input power to the titanium model and prototype was extremely close 
to that which the parameters required, which implied a negligible error 

Miller (21) discussed a means of estimating thermocouple loss. The 
thermocouple leads were considered to be infinitely long pm fins with a 
known root temperature, and with radiation heat transfer from the fm 
surface. The differential equation that describes this situation is* 
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( 76 ) 


This equation may he integrated once with respect to x to yield the 
heat transfer rate at the root of the fm, the desired quantity. 


= -kA (~| = [*- n 2 R 3 kea (T 5 - T 5 ) 

nldx / „ 15 o ® 

\ / x = 0 L J 


The wire radii were known from tabulated values for 30 gage wire. 
Thermal conductivities of 220 BTU/hr - ft - °F and 14 BTU/hr - ft - °F 
for copper and constantan, respectively, were used in calculations. The 
surface emittance values were estimated for the polished wires as 0.05 for 
copper and 0.10 for constantan. Measured values were used for T q and T^ 

T7as assumed to be 150°B., the boiling temperature of liquid nitrogen. 

The total thermocouple loss was less than one per cent. 


Losses through the nylon support lines were considered to be negli- 


gible . 

The problem of pover lead loss (or gam) was eliminated by scaling 
the leads using the technique m Appendix A. A dramatic temperature drop 
from the center to the ends of the rods was observed when oversized leads 
were used. The reverse trend was present when undersized leads were used 
and acted as heat sources for the rods. 

In the derivation of the similarity parameters, it was assumed that 
model and prototype were constructed of the same material, which implies 
that thermal conductivity and specific heat were equal m model and pro- 
totype. Since the percentage of alloying elements for a particular material 
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xs not the same for each batch from a particular manufacturer or for 
materials with the same alloy specifications from different manufacturers, 
there was a possibility that model and prototype had slightly different 
thermal properties. 

Miller (21) assumed that material properties varied with temperature 
according to the relationships* 

k = k T a (78) 

o 


c 

p 


c 

po 


(79) 


The effects of the variation* from lot-to-lot, of the conductivity and 

specific heat on the value of the modeling parameters was illustrated by 

assuming that the exponents in his power la w equations were the same for 

both model and prototype which determined the effect of the variations of 

k and c , Miller found that a three per cent variation m k could 
o po r o 

cause a four percent variation in the input energy parameter, a three 
per cent variation xn the temperature parameter. Comparison of resistance 
per foot data for models and prototypes indicated that a similar trend 
existed xn this investigation. 
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CONCLUSIONS 


The temperature-time data m Figure 8 and 9 verify the applicability 
of thermal modeling to a transient system nth internal generation. The 
objective of the investigation ras satisfied by experimental verification 
of the analytically derived similarity parameters. The development of a 
power lead scaling technique and presentation of error analysis techniques 
should be useful m future thermal modeling endeavors. 
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APPENDIX A 


TECHNIQUE OF SIZING POWER LEADS 



TECHNIQUE OF SIZING POWER LEADS 


The problem of heat loss due to conduction at the ends of the 
prototypes and models proved to be a problem of significant magnitude 
during the experimental phase of the investigation. After many test 
runs, it was determined that the size of the power leads connected at 
each end of a model or prototype was the chief source of this problem. 
If the leads were much larger m diameter than the test rods, they 
served as heat sinks at the rod ends. However, under sizing the leads 
tended to reverse the situation, and a heat source was created at the 
rod ends. It was thus necessary to develop a technique of scaling the 
power leads for each test rod m such a way that the leads would not 
serve as a heat sink or heat source. 


The power supplied to a test rod can be expressed as: 

P = i 2 r (1) 

Since the same current passes through the leads as through the rod, 
the following equation must hold: 


Zl Zl 

p r _r r 


( 2 ) 


The above resistances have the unit ohms per foot of length. 

The power supplied to each test rod is dissipated by radiation at 
the surface of the rod. It is evident, then, that the watts of power 
supplied can be expressed m terms of heat radiated. The ratio of the 
heat radiated from the lead and rod is given by: 
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( 3 ) 



Vl t l a 


a R £ R T R A 


sL 


sh 


The emxssivity terms were made equal by painting the leads black. 

Thxs had the effect of keeping the lead cooler, since it closely 
approached being a perfect radiator. Conversely, polishing the leads 
tended to keep them warmer. For steady state temperature between the 
leads and rod, equation (3) becomes: 




(4) 


Since heat radiated and power 
(.293 watts/BTU/hr) , equations (2) 


are related by a constant conversion, 
and (4) can be related to give 


A sL = 5k 
A sR \ 


(5) 


Equation (5) reduces to. 
°R \ 


( 6 ) 


The resistance per foot of length of the rod can be found by running 
a known current through the lead and measuring the voltage drop with a 
voltmeter* The application of Ohm T s law will then give the test rod 
resistance The diameter of the rod was 

be reduced to a constant C 5 and equation (6) becomes - 

D l = (7) 


already known, so, D^/R^ may 
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Equation (7) xs obvxously lxnear. 

To utxlxze equatxon (7), a plot of lead dxameter xn xnches vs re- 
sxstance in ohms per foot of length was made by usxng a table of the 
propertxes of copper wxre. By pxckxng several values of xn equatxon 
(7) , xt was possible to draw a straight line through the above mentioned 
curve, where the origin of this line began at = 0 and = 0. This 
starting point was found by simply substituting - 0 into equation (7), 
which resulted m D =0. The intersection of the curve and the line 

1j 

given by equation (7) gives the proper diameter to use for each test rod. 
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APPENDIX B 


Conversion of Millivolt Readings 
to Temperature Readings 



Conversion of Millivolt Readings 
to Temperature Readings 

Some explanation is m order concerning the procedure followed m 
converting the millivolt reading from the recorders to temperature readings. 
Holman (24) stated that it is common to express thermoelectric emf m terms 
of the potential generated with a reference junction at 32 °F. Standard 
thermocouple tables have been prepared on this basis showing the output 
characteristics of many common thermocouple combinations. Such a table 
was used for the copper-constantan thermocouples utilized in this experiment. 

A Leeds and Northrup 8690 potentiometer was connected directly to 
each strip-chart recorder, and the slidewire scale setter was turned on 
the potentiometer until the pointer on the strip-chart recorder was at 
the zero point on the recorder paper. The strip-chart paper was divided 
into one-hundred intervals. When the pointer was above the zero point 
on the recorder paper, the millivolt reading on the potentiometer was 
recorded. This procedure was followed m intervals of ten across the 
recorder paper until the one-hundred intervals had been crossed. Then 
the copper-constantan thermocouple table referenced at 32°F was used to 
match the data m the table with the data recorded. It was found that 
the recorders had a tolerance of + 1.0 degree Fahrenheit. From the re- 
corded data a plot of percent of millivolt range across recorder paper 
vs. temperature was made as shorn in Figure B-l. This curve was used 
to determine model and prototype temperature from the strip chart recorders. 
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